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An inhomogeneous backflow transformation for many-particle wave functions is presented and 
applied to electrons in atoms, molecules, and solids. We report variational and diffusion quantum 
Monte Carlo (VMC and DMC) energies for various systems and study the computational cost of 
using backflow wave functions. We find that inhomogeneous backflow transformations can provide 
a substantial increase in the amount of correlation energy retrieved within VMC and DMC cal- 
culations. The backflow transformations significantly improve the wave functions and their nodal 
surfaces. 
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I. INTRODUCTION 

The fermion sign problem continues to preclude the apphcation of in principle exact quantum Monte Carlo (QMC) 
methods to large systems, and so approximate QMC methods must be used instead. Probably the most widely-used 
of these is the stable and efficient diffusion quantum Monte Carlo (DMC) algorithmj^ll^l in which the fermion sign 
problem is sidestepped through the introduction of the fixed-node approximation.^ DMC can provide highly accurate 
energies for assemblies of quantum particles, but the fixed-node approximation is uncontrolled and its accuracy is 
often difficult to assess. 

The fixed-node approximation'^ involves constraining the nodal surface of the wave function to equal that of an 
approximate "trial" or "guiding" wave function. The fixed-node DMC energy is higher than the ground-state energy, 
becoming equal in the limit that the fixed nodal surface is exact. The dependence of the DMC energy on the quality 
of the trial wave function is often significant in practice. It would therefore be very useful to be able to construct trial 
wave functions with better nodal surfaces to reduce the effect of the fixed-node approximation. 

Efforts to construct wave functions with accurate nodal surfaces have continued since the introduction of the fixed- 
node approximation. Single-determinant wave functions often provide good nodal surfaces for closed-shell systems, and 
multideterminant wave functions can do so for small open-shell systems, although the required number of determinants 
becomes excessive for large systems. Compact pairing wave functions consisting of an antisymmetrized product of 
two-electron "geminals"l3 were introduced long agcP^and have recently been used in QMC calculations for atoms and 
molecules .'^ Triplet-pairing Pfaffian wave functions were first used in QMC calculations for liquid ^He by Bouchaud 
and Lhuillier,^ and recently this approach has been extended by Bajdich et a/.,^*^ who considered atomic and molecular 
systems in which both parallel- and antiparallel-spin electrons are paired. 

Another approach for improving upon a single determinant of one-electron orbitals is to introduce parameters 
which allow the orbitals to depend on the positions of the other electrons. Such a route was followed by Wigner and 
Seitz,'li^ who considered wave functions in which the orbitals of the up-spin electrons depend on the positions of the 
down-spin electrons, and vice versa. This idea surfaced again much later in connection with the quantum-mechanical 
description of "backflow." Classical backflow is related to the flow of a fluid around a large impurity. Its quantum 
analog was discussed by FeynmarP^ and Feynman and CoherP^ in the contexts of excitations in ''He and the effective 
mass of a ^He impurity in liquid ^He. They argued that the energy would be lowered if the ^He atoms executed 
a flow pattern around the moving '^He impurity which prevented the atoms overlapping significantly. This effect 
was shown to correspond to the requirement that the local current of particles is conserved. They recognized that, 
without backflow, the effective mass of the '^He impurity would equal the bare mass, and incorporating backflow led 
to a substantial increase in the effective mass. It turns out that the mathematical form obtained by incorporating 
backflow into a single-determinant wave function is related to the wave functions considered by Wigner and Seitz.'ii' 

In later studies, wave functions including Jastrow factors and backflow-like correlat ions we re used to study a '^He 
impurity in liquid '*He and liquid '^He within a Fermi-hypernetted chain approximation.'^lMiS] Backflow was first used 
in QMC calculations by Lee et al.^ who calculated the total energy of liquid '^He. QMC calculations for electrons 
using Slater- Jastrow wave functions with backflow correlations were first performed by Kwon et al.^^ for the two- 
dimensional homogeneous electron gas (HEG) , and latei'ISi for the three-dimensional HEG (see also the paper by Zong 
et al^. QMC calculations using Slater-Jastrow wave functions with backflow correlations have also been reported 
for solid and liquid hydrogeUjl^IISSI which were the first such applications to inhomogeneous electron systems. 

While Jastrow factors keep electrons away from one another and greatly improve wave functions in general, they 
do not alter nodal surfaces. Holzmann et aZ.l^have argued that backflow and three-body Jastrow correlations arise 
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as the next-order improvements to the standard Slater- Jastrow wave function, which consists of a Slater determinant 
multiplied by a two-body Jastrow factor. The importance of backfiow correlations within DMC calculations is that 
they alter the nodal surface and can therefore be used to reduce the fixed-node error. 

In this paper we introduce parameterized inhomogeneous backfiow transformations, and apply them to atoms, 
molecules, and extended systems. The rest of this paper is structured as follows: a general description of the Slater- 
Jastrow and backfiow wave functions is given in Sec. Ill an explicit form for the backfiow displacement field is developed 
in S ec.|III[ an extensive set of results is given in Sec. Wjand discussed in Sec.|V] and our conclusions are summarized in 
Sec. |VI[ Iniportant technical information about the calculations, including the constraints on the backfiow parameters, 
has been gathered in the appendices. Hartree atomic units {h ^ \e\ = me 47reo = 1) are used throughout. 



II. SLATER-JASTROW AND SLATER-JASTROW-BACKFLOW WAVE FUNCTIONS 



The Slater- Jastrow (SJ) wave function can be written as 

vI'|J(R) = e'^(^)vl/s(R) , (1) 

where R denotes the set of electron coordinates {r^}, e'^'^^^ is the Jastrow correlation factor, and the Slater part 
^'s(R) consists of a determinant or sum of determinants, defining the nodes of ^|.''(R). 

Backfiow (BF) correlations are introduced by substituting a set of collective coordinates X for the coordinates R 
in the Slater determinants, so that 

*|P(R) = e^(^%s(X) , (2) 
where each of the new coordinates is given bjff^l^Il 

X, =r, + ^,(R), (3) 
where is the backfiow displacement of particle i, which depends on the configuration of the whole system. 



III. INHOMOGENEOUS BACKFLOW TRANSFORMATIONS 

The form of the backfiow displacement in homogeneous systems has been taken a^^^'^^'^H 

Vij^ij J (4) 



We 



where is the number of electrons and rjij = ri(rij) is a function of the interparticle distance r.^. Eq. Q is the 
most general isotropic two-electron coordinate transformation for a homogeneous system. A single electron i perceives 
space to be isotropic, but when another electron j is introduced, the electron-electron (e-e) vector Vij becomes an 
inequivalent direction. The e-e backfiow displacement is taken to be along this direction, as there is no reason why a 
displacement in a specific perpendicular direction should occur. 

In a system with nuclei a new set of directions is introduced, the electron-nucleus (e-n) vectors {fi/}, and one is 
led to introduce an c-n contribution to ^j, of the form 



fen ^ 
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A*i/ri7 , (5) 



where jiu — f^irn) and Nn is the number of nuclei. 

We also introduce an electron-electron-nucleus (e-e-n) term to describe two-electron backfiow displacements in the 
presence of a nearby nucleus, 

We N„ 

where <^{^ — {rn ,rjj ,rij) and 6^^ — Q^{rii,rjj,rij). Note that the vector ^{^Vij + ^r^/ is capable of spanning 
the plane defined by r,, r^, and r/, without the need to introduce a component along the direction of Vjj. The total 
backfiow displacement is the sum of these three components, = ^""^ -I- ^™ -I- 
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At large distances 'ri{rij) is expected to decay as r~-^ in three dimensionJl^l and t^^'^'^ in two dimensions However, 
for computational efficiency and for compatibility with periodic boundary conditions it is better to cut off the 77 
function and the other backfiow functions smoothly at some radius. We use a simple cutoff function, 



L — r 



f{r-L)=^-j- j H{L~r) , (7) 

where r is to be substituted by an e-e or e-n distance as appropriate, L is the cutoff length, C is the truncation 
order^ and H denotes the Heaviside function. The advantages of this cutoff function are, firstly, that its value can 
be computed rapidly and, secondly, that it has considerable flexibility because one can choose the value of C and use 
L as an optimizable parameter. 

Rational^ and Gaussian'^^ forms for homogeneous backflow functions have been used in previous work. However, 
we have chosen to use natural power expansions because of the excellent results we have obtained with such expansions 
for our Jastrow factor,!^ and the lack of a priori knowledge of more specific parameterizations for the inhomogeneous 
functions. It is estimated that numerical errors in the evaluation of natural polynomials become significant beyond 
order about 20 when using double-precision arithmetic and, although one can go to substantially larger orders using 
Chebyshev polynomials, we have not found this to be an issue in our work. 

We have used the following polynomial expansions for r/ij, fx^j, ^j^, and 6^^, 



k=0 

fJ-ri = f{ni;L^j)J2dkjr'li, (9) 

fc=0 

A^cn,! Afon.f ^cc,/ 

i.f = f{ru;Uj)f{r,r,L^j) E E E VUmjr'yjr'.ir'^, (10) 

fe=o ;=o m=0 

N^r^J N,n,I A^cc.J 

ef = f{ru;Uj)f{r,i;L^j) E E E (^kimjr'^jr'.irT;, (11) 



fe=0 1=0 m=0 



where N^i, N^j, N^nj, and Nccj are the expansion orders, i,,, L^j, and L^j, are cutoff lengths, and {cfc}, {dk,i}, 
{fkimj}, and {9kim,i} are the optimizable parameters. We allow the parameters in rj, <i>^, and 6^ to depend on the 
spins of the electron pairs, and those in fi to be spin dependent; for simplicity we have omitted such dependencies in 
the description of the functional forms above. In periodic systems, we constrain and to be smaller than the 
Wigner-Seitz radius iws of the simulation cell and L$ to be smaller than for computational efficiency. 



IV. RESULTS 



In this section we present variational quantum Monte Carlo (VMC) and DMC results obtained with our imple- 
mentation of backffow transformations. The CASINO cod4^ has been used for all of our QMC calculations. Our 
DMC algorithm is essentially as described in Ref. [25] All DMC energies reported here have been extrapolated to zero 
time step. We have optimized the parameters in our wave functions by minimizing the unreweighted va rianc e of the 
energy^ using a scheme which facilitates the optimization of parameters that modify the nodal surface .12^1221 

We have used the Jastrow correlation factor of Drummond et al.^ In our all-electron (AE) calculations, with the 
exception of those for the HEG, the orbitals were obtained from Hartree-Fock (HF) calculations using large Gaussian 
basis sets and the CRYSTAl98 code,'^ and the cusp-correction algorithm of Ref. [31] was applied to each orbital 
at each nucleus. In our pseudopotential (PP) calculations we used the Dirac-Fock Average Relativistic Effective 
pseudopotentials of Refs. 1^ and [551 the nonlocal energies being calculated within the locality approximation.^'* The 
one-electron orbitals were obtained from the plane- wave PP CASTEP cod#^ using the Perdew-Burke-Ernzerhof (PBE) 
generalized-gradient-approximatioiPSl exchange-correlation functional. The orbitals were re-expanded in terms of 
"blip" functions,'^' making the QMC calculations much more efficient. 

We have reported the variance of the total local energy for our VMC calculations, = {H'^) — (H)'^,^^ while the 
reported mean energies are either total, per electron, or per primitive cell, as we have found appropriate in each case. 
We have estimated the amount of correlation energy retrieved in our calculations by comparing our energies with 
"exact" reference data, where available. In the case of the PP carbon atom and PP carbon dimer we have used the 
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estimates of the PP valence correlation energy of Ref. [351 assuming an error bar of 0.004 a.u. as suggested by the 
author. In the HEG we have used our BF-DMC energies as if they were "exact" , and in PP carbon diamond we have 
not estimated the amount of correlation energy retrieved. 

A. Homogeneous electron gas 

We studied three-dimensional, unpolarized HEGs consisting of 54 electrons in a simple cubic simulation cell subject 
to periodic boundary conditions. As well as the densities of Ts = 1, 5, 10, and 20 studied by Kwon et al^ and 
Holzmann et aZ.'^ using backflow wave functions, for completeness we studied two additional densities, = 0.5 and 
2. Holzmann et al. used an analytical backflow form containing no variable parameters in addition to a Gaussian form 
with variable parameters. In each case we compare our result with the corresponding lowest-energy backflow result 
from table II of Ref. [HI 

We included a plane- wave term in our Jastrow factor [Eq. (28) of Ref. [24], which we found to improve the variational 
energies at all densities. We also studied the effect of including a symmetric three-electron Jastrow term, W , of the 
type used in Ref. UHl with 

= X! X! X! (^^i'^u) ■ iwikVik) , (12) 

where Wij is a function of the distance between electrons i and j, which we parameterized as 

Wij = firij-, Lu,) ^ eir\j , (13) 

where Nyj is the order of the expansion, the {e;} are expansion parameters, and / is the cutoff function of Eq. ([t]). We 
decided to include a W term for all densities at the Slater- Jastrow level, while we used it in conjuction with backflow 
only for the three lowest densities, where its effect on the SJ energy was found to be statistically signiflcant. We refer 
to the SJ and BF wave functions with a three-electron Jastrow term as SJ3 and BF3, respectively. The backflow 
parameters were allowed to depend on the spins of the electron pairs, while the parameters in the three-electron 
Jastrow factors were constrained to be independent of spin, as this gave slightly better results. The expansion orders 
and iV^, were set to 8 for all densities. The cutoff lengths L,, and were optimized, but at all densities they 
adjusted themselves to the maximum allowed value (the Wigner-Seitz radius). The energies and variances obtained 
are given in Table |T] and the energies are illustrated in Fig. [T] which gives the percentage of the correlation energy 
retrieved at different levels as well as the SJ and BF energies of Ref. [ST] The introduction of backflow increases the 
kinetic energy, but decreases the potential energy by a larger amount. Our SJ-DMC energies are in good agreement 
with those of Holzmann et al., which of course they should be, because the SJ trial wave functions have identical nodal 
surfaces. Our SJ-VMC calculations retrieve a higher percentage of the correlation energy than those of Holzmann 
et al., and we believe this is mainly due to the plane- wave term in our Jastrow factor. Our BF-VMC calculations 
consistently retrieve 99.5% of the correlation energy throughout the density range considered, while those of Holzmann 
et al. drop below 99% for > 5. Our BF-DMC energies are within error bars of those of Holzmann et al. In agreement 
with the work of Refs. [2T|and[iP, we found that backflow gives a larger energy reduction at the VMC level than the 
three-body Jastrow term W at all densities, although W becomes more important at large r^. 

The variances of the VMC energies reported in Table |l] are illustrated in Fig. [2] The lines on the log-log plot 
corresponding to our SJ and BF variances are almost parallel, indicating an almost constant ratio of the SJ to BF 
variances of about 4. The variances of Holzmann et al. are systematically higher than ours for comparable calculations, 
and at — 20 our SJ variance is lower than their BF variance.!^ 

The optimized homogeneous backflow displacement rjirii)rij is plotted in Figs. [3]and|4] and the optimized three- 
body function Wij is shown in Fig. [s] Holzmann et al^^ and Kwon et al^^ used identical rj functions for parallel 
and antiparallel spin pairs, whereas we have allowed them to differ. At each density, the maximum value of 'ri{rij)rij 
for antiparallel spins is over twice as large as that for parallel spins, and occurs at smaller electron separations. The 
backflow displacements for antiparallel spins are generally larger than for parallel spins, and hence antiparallel-spin 
backflow is much more important than parallel-spin backflow. Our antiparallel-spin 77 function is similar to the spin- 
independent rj function of Kwon et al.^^ except that we do not find an attractive tail at — 20. Note that, to 
obey the cusp conditions, we constrain the parallel-spin rj{ri.j) function to have zero derivative at Vij = 0, while the 
antiparallel-spin 77 function may have a nonzero derivative: see Appendix ] A 1| This accounts for the differences in the 
behavior of the parallel- and antiparallel-spin 77 functions at small Vij which are visible in Figs. [3] and [4] 
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rs Wfn. Ey (a.u./elec) 


o-y (a.u.) 


CEv (%) (a.u./elec) 


CEd (%) 


0.5 HF 


3.2659(7) 


76(1) 


0(3) 


- 




SJ 


3.2236(2) 


3.34(3) 


94.5(5) 


3.22245(9) 


97.0(3) 


SJ3 


3.2233(2) 


3.4(2) 


95.1(5) 






Ht 


3.22132(7) 


0.76(1) 


99.5(2) 


3.22112(4) 


100.0(2) 


1.0 HF 


0.5689(4) 


19.1(4) 


0(2) 


- 


- 


SJ 


0.53211(7) 


0.719(7) 


94.3(3) 


0.53089(9) 


97.5(4) 


SJ3 


0.53175(7) 


0.80(6) 


95.3(3) 








0.o3U0y(3j 


0.163(2) 


yy.o^z ) 


0.02989(4) 


100.0(2) 


2.0 HF 


0.0186(2) 


4.9(1) 


0(1) 


- 


- 


SJ 


-0.01246(3) 


0.147(2) 


95.4(1) 


-0.01311(2) 


97.4(1) 


SJ3 


-0.01252(3) 


0.138(2) 


95.6(1) 






at 


—0.01382(2) 


0.0342(6) 


c\c\ T^a /'t\ 

99.56(7) 


C\ C\'\ C\Ci O ( C\\ 

—0.013966(9) 


1 r\c\ c\r\ ( ci\ 

100.00(6) 


5.0 HF 


-0.05625(7) 


0.76(1) 


0.0(6) 


- 


- 


SJ 


-0.07815(1) 


0.0149(2) 


96.09(7) 


- 




SJ3 


-0.078284(9) 


0.0129(3) 


96.70(6) 


-0.078649(7) 


no on / \ 

98.30(5) 


±3t 3 


—0.078961(5) 


0.00317(6) 


c\c\ t/ o \ 

99.67(3) 


—0.079036(3) 


1 nn nn /" o \ 

100.00(3) 


10.0 HF 


-0.03884(4) 


0.194(4) 


0.0(5) 


- 


- 


SJ 


-0.053927(4) 


0.00236(2) 


96.69(4) 


- 




SJ3 


-0.054042(4) 


0.00179(3) 


97.43(4) 


-0.054255(3) 


98.80(4) 


BF3 


-0.054389(2) 


0.00055(1) 


99.65(2) 


-0.054443(2) 


1 nn nn /nN 

100.00(2) 


20.0 HF 


-0.02205(2) 


0.0477(9) 


0.0(4) 






SJ 


-0.031767(2) 


0.000377(4) 97.20(4) 






SJ3 


-0.031858(1) 


0.000237(2) 


98.11(3) 


-0.031973(3) 


99.26(5) 


BF3 


-0.0319984(8) 


0.000091(1) 99.51(3) 


-0.032047(2) 


100.00(3) 



TABLE I: Energies and variances for three-dimensional, unpolarized HEGs consisting of 54 electrons in a simple cubic simulation 
cell. _Ev and Et, refer to VMC and DMC energies, respectively; CEy and CEd are the percentages of the correlation energies 
retrieved at the VMC and DMC levels, respectively, and a\ is the VMC variance. 



The magnitude of our optimized three-electron Jastrow factor, represented in Fig. |5] increases monotonically with 
Ts, and the maximum of 6(rij'Wij)'^ is at about Tijjrs = 0.4 for all densities. This is in contrast with the behavior of the 
three-electron Jastrow factor of Kwon tt al. (see Fig. 1 of Ref . . which changes sign at = 1 (our parametrization 
is not allowed to do so) and breaks its monotonicity with rg at = 20. Kwon et al. find that the maximum of the 
plotted function is located at about rij/rg — 1. 

B. Lithium atom and dimer 

1. AE lithium atom 

Our results for the ground state of the AE lithium atom are given in Table |ll] The SJ wave function gives a 
reasonably good VMC energy. Our backflow function consists of a spin-pair-dependent e-e-n term with Nen = N^e — 3; 
this produces a BF-VMC energy that is within statistical error bars of the exact value. Note that the BF-VMC, SJ- 
DMC, and BF-DMC energies are within statistical error bars of each other and are very close to the exact value. The 
excellent performance of the BF-VMC calculation is particularly noteworthy. The single-determinant nodal surface 
of the ground state of lithium is certainly extremely accurate and may even be exact, although some contrary 
evidence has been cited. '^^ It is therefore unlikely that backflow could improve upon the SJ-DMC energy, and indeed 
it leaves it essentially unchanged. 
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FIG. 1: (Color online) Percentages of the correlation energy recovered at the (solid circles from top to bottom) BF-DMC, BF- 
VMC, SJ-DMC, and SJ-VMC levels as a function of the density parameter (see Table[l|. Zero correlation energy corresponds 
to HF-VMC and 100% to BF-DMC. The hollow circles are the best BF-DMC, BF-VMC, SJ-DMC, and SJ-VMC energies of 
Holzmann et al.^^'m the same order. The statistical error bars on the QMC data are smaller than the symbols except where 
error bars are shown. 




FIG. 2: (Color online) VMC variances achieved at the SJ level (top solid line) and BF level (bottom solid line) as a function of 
the density parameter (see Table The hollow circles are the best SJ-VMC and BF-VMC variances of Holzmann et al.^^ 
in the same order. The statistical error bars on the QMC data are smaller than the symbols except where error bars are shown. 
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FIG. 3: (Color online) Antiparallel-spin homogeneous backflow displacement rj{rij)rij for the HEG at the different densities 
studied. For the three highest densities, the curves correspond to BF wave functions, while the others are for BF3 wave 
functions. 
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FIG. 4: (Color online) Parallel-spin homogeneous backflow displacement r]{rij)rij for the HEG at the different densities studied. 
For the three highest densities, the curves correspond to BF wave functions, while the others are for the BF3 wave function. 



2. AE lithium dimer 



We studied the ground state of the AE Li2 dimer at the experimental bond length of 5.051 a.u.l^ We tested several 
different backflow functions, obtaining the results given in Table [ITT] The use of homogeneous backflow retrieves only 



8 



1.0 



0.8 



0.6 



^ 0.4 



0.0 



0.2 h / 

7 •■ 



N \ 



r =0.5 

s 

-- r=\ 

s 

.-. r=2 
.... r=5 

/=10 

.-. r=20 

s 



0.0 



1.0 



2.0 



r.. I r 



3.0 



FIG. 5: (Color online) Three-body contribution to the Jastrow function for the HEG due to three electrons at the vertices of 
an equilateral triangle of side Vij, at the different densities studied. For the three highest densities, the curves correspond to 
SJ3 wave functions, while the other three are for BF3 wave functions. 



Method Wfn. Np E (a.u.) (a.u.) % corr. en. 



HF 






-7.43273 


0.0 


Exact 






-7.47806 


100.0 


VMC 


SJ 





-7.47648(3) 0.00385(2) 


96.52(8) 




BE 


114 


-7.47801(3) 0.00241(1) 


99.89(6) 


DMC 


SJ 





-7.47803(8) 


99.9(2) 




BE 


114 


-7.47802(6) 


99.9(1) 



TABLE II: Slater- Jastrow and backflow results for the AE lithium atom. The number of free backflow parameters, excluding 
cutoff lengths, is Np. The Hartree-Eock (HE) and exact energies were taken from Refs. BOl and 1411 



an additional 0.7% of the correlation energy. A plot of the VMC energy as a function of the number of parameters 
is displayed in Fig. [6] which shows that the reduction in VMC energy is very small beyond about 150 parameters. 
Whereas backflow gave 99.89(6)% of the correlation energy at the VMC level for the lithium atom, for the dimer our 
best backflow transformation retrieves only 87.79(8)%. At the DMC level the improvement is small: using a SJ wave 
function we obtain 96.2(3)% of the correlation energy while with the backflow wave function this improves slightly [to 
97.1(3)%]. Considerably better DMC results for Li2 have been obtained using multideterminant (MD) wave functions. 
Bressanini et a/. ''^ obtained a DMC energy of —14.9923(2) with one configuration state function (CSF), while their 
best result was -14.9952(1) with 4 CSFs. 

The computed binding energies of Li2 are presented in Table [rv] The BF-VMC, SJ-DMC, and BF-DMC energies 
of the AE lithium atom are within error bars of the exact energy, and therefore the error in the binding energy arises 
solely from the Li2 energy. Backflow improves the VMC and DMC binding energies of Li2 a little, but it is still 
somewhat short of the exact value. The single-determinant nodal surface of Li2 is quite poor, and backflow is not 
very effective at improving it. Combining MD wave functions with backflow might yield signiflcant improvements in 
this case. 
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Method Wfn. 






iVcn 


iVcc 


Np 


E (a.u.) 


(a.u.) 


% corr. en. 


HF 




_ 




_ 


_ 


_ 


-14.871545 


_ 


0.0 


Exact 


- 












— 14.9954 




100.0 


VMC 


SJ 


- 


- 


- 


- 





-14.9751(1) 


0.0165(1) 


83.6(1) 




BF 





6 








14 


-14.9755(1) 


0.01607(9) 


83.9(1) 




BF 


8 











17 


-14.9760(1) 


0.01590(7) 


84.3(1) 




BF 








2 


2 


16 


-14.9768(1) 


0.01424(7) 


84.9(1) 




BF 








2 


4 


44 


-14.9782(1) 


0.01273(7) 


86.15(9) 




BF 








2 


6 


72 


-14.9789(1) 


0.0125(1) 


86.65(9) 




BF 








3 


4 


156 


— 14.9797(1) 


0.01102(5) 


87.33(8) 




BF 


n 

u 


n 

u 




A 


ouo 










BF 








4 


3 


230 


-14.9803(1) 


0.01038(4) 


87.79(8) 


DMC 


SJ 













-14.9907(4) 




96.2(3) 




BF 








4 


3 


230 


-14.9918(4) 




97.1(3) 



TABLE III: Slater-Jastrow and backflow results for the AE Li2 molecule. The different backflow forms have been put in order 
of decreasing energy. The number of free backflow parameters, excluding cutoff lengths, is A'p. The Hartree-Fock (HF) and 
exact energies were taken from Ref. 1441 
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FIG. 6: The VMC energy of AE Li2 versus the total number of parameters in the backflow functions. 



C. Carbon atom, carbon dimer, and diamond 

1. AE carbon atom 



The ground state of the AE carbon atom is a good example of a system where single-determinant wave functions 
result in large fixed-node errors: see Table [V] In this case, we have tested several combinations of terms, expansion 
orders, and constraints to explore the possibilities of backflow transformations. The VMC data in Table |V) and 
additional data, are plotted in Fig. [7] where the performance of the different backflow functions used can be compared 
conveniently. Using only homogeneous backflow (flrst BF-VMC results in Table lv| gives a very small reduction 
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Method Wfn. Et (a.u. 



HF 




0.0061 


Exact 




O.Ooao 


VMC 


SJ 


0.0221(1) 




BF 


0.0243(1) 


DMC 


SJ 


0.0346(4) 




BF 


0.0358(4) 



TABLE IV: Slat er-Jastrow and backflow binding energies for the AE Li2 molecule, computed using the best results from 
Tables In] and III The Hartree-Fock (HF) and exact energies were taken from Refs. 1401 1411 and 1441 and references therein. 



in energy. It seems that inhomogeneous systems require inhomogeneous backflow to produce good wave functions, 
and the e-e-n term is particularly successful in providing this. To evaluate the relative importance of the two e-e-n 
functions and 6^^, we performed calculations constraining the parameters in one of them to be zero. The results 
are also given in Table [v| In this case Qj^, which contributes to the backflow displacement in the direction of ru, 
is slightly more important than . Applying both terms gives better results than using only one of them, as we 
expected. We also tested the effect of constraining the backflow displacement to be irrotational, which was suggested 
in Ref . 1211 The application of this constraint , which is explained in Appendix |A 2[ approximately halves the number 
of parameters in the backflow functions, but it gives very poor results for the carbon atom. 



Method Wfn. 








S 


I 




E (a.u.) 


(a.u.) 


% corr. en. 


HF 
















-37.688619 




0.0 


Exact 
















-37.8450 




100.0 


VMC 


SJ 















-37.8064(3) 


0.193(2) 


75.3(2) 




BF 


8 








T 




17 


-37.8089(3) 


0.194(2) 


76.9(2) 




BF 





6 





T 




10 


-37.8089(3) 


0.184(1) 


76.9(2) 




BF 








2 


F 


F 


10 


-37.8119(3) 


0.1685(8) 


78.8(2) 




BF 


8 


6 





T 




27 


-37.8126(3) 


0.178(2) 


79.3(2) 




BF 








4 


T 


T 


35 


-37.8140(3) 


0.171(2) 


80.2(2) 




BFt 








3 


T 


F 


56 


-37.8155(3) 


0.1578(9) 


81.1(2) 




BF 








5 


T 


T 


121 


-37.8177(3) 


0.153(2) 


82.5(2) 




BF 








2 


T 


F 


16 


-37.8180(3) 


0.159(4) 


82.7(2) 




BE* 








3 


T 


F 


58 


-37.8198(3) 


0.144(2) 


83.7(2) 




BF 








3 


F 


F 


60 


-37.8225(3) 


0.135(2) 


85.5(2) 




BF 








4 


F 


F 


158 


-37.8239(3) 


0.119(6) 


86.5(2) 




BF 








3 


T 


F 


114 


-37.8246(3) 


0.127(2) 


87.0(2) 




BF 





6 


3 


T 


F 


124 


-37.8252(3) 


0.122(1) 


87.3(2) 




BF 








4 


T 


F 


308 


-37.8259(3) 


0.109(1) 


87.8(2) 


DMC 


SJ 















-37.8297(2) 




90.2(1) 




BF 





6 


3 


T 


F 


124 


-37.8324(1) 




92.0(1) 



TABLE V: Slater-Jastrow and backflow results for the AE carbon atom. The different backflow forms have been put in order 
of decreasing energy. Key: A^^ = Ncn ~ Nac', S indicates whether the parameters are spin and spin-pair dependent (T) or not 
(F); I indicates whether the constraints for irrotational backflow have been applied (T) or not (F); Np is the number of free 
backflow parameters, excluding cutoff lengths. Where the f symbol is used, we constrained 6kim,i = 0; where the ^ symbol 
appears, we constrained ifikimj ~ 0. The Hartree-Fock (HF) and exact energies were taken from Refs. BDl and 1411 



The most satisfactory backflow forms reduce the difference between the VMC and exact energies by a factor of 
about 2. The further energy reduction from using DMC is quite small, and our BF-DMC calculation gave an energy of 
—37.8324(1) a.u., which corresponds to 92.0(1)% of the total correlation energy. This suggests that, although backflow 
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FIG. 7: (Color online) The VMC energy versus its variance for AE carbon (see Table \ V\ . The open symbols denote that 
the backflow parameters are indepen dent of spin, while the filled symbols denoted spin-dependent parameters. The exact 
nonrelativistic energy is —37.8450 a.u.'^^'^ 



improves significantly upon the single-determinant nodal surface of the carbon atom, it misses some important features 
of the exact nodal surface. It is well-known that the single-determinant nodal surface of the carbon atom can be 
substantially improved by using MD trial wave functions. Barnett et ali^ used an MD trial wave function consisting 
of 14 CSFs, and obtained a DMC energy of —37.8420(3) a.u., which corresponds to 98.1(2)% of the correlation energy. 
Glauser et aim showed that the configuration space of a single-determinant of HF orbitals for the ground state 
carbon atom is divided into four nodal pockets,'^ but more accurate wave functions indicate that the exact wave 
function has two nodal pockets. It appears that backfiow transformations are unable to correct this defect in the 
single-determinant nodal surface. 



2. PP carbon atom 



We have also studied how backfiow performs in systems where PPs are used. Our results for a PP carbon atom are 
given in Table |VI[ The reduction in the VMC energy obtained with backfiow is much smaller than for the AE carbon 
atom, but the corresponding energy reduction within DMC of 0.0039(1) a.u. is somewhat larger than the AE one of 
0.0027(2) a.u. A peculiarity of this case is that the reduction in the DMC energy resulting from the use of backfiow 
is 71% of the reduction in the VMC energy, which is the largest such percentage amongst the calculations described 
here. 



3. PP carbon dimer 

For the PP carbon dimer we used the experimental bond length of 2.3622 a.u.,!^ obtaining the results given in 
Table |VII[ The carbon dimer is another example of a system in which MD effects are known to be substantial. 
Backflow results in larger energy reductions per atom than for the isolated atom at both the VMC and DMC levels. 
The computed binding energies of C2 are presented in Table |VIII[ The use of backflow slightly improves the binding 
energy of the dimer. 
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Method Wfn. Np E (a.u.) (a.u.) % corr. en. 



HF 






-5.31663 




0.0 


Exact 






K A on/ A ^ 






VMC 


SJ 





-5.4007(1) 


0.0582(4) 


81(3) 




BF 


218 


-5.4061(1) 


0.0502(6) 


87(3) 


DMC 


SJ 





-5.40886(7) 




89(4) 




BF 


218 


-5.41273(9) 




93(4) 



TABLE VI: Slater-Jastrow and backflow results for the PP carbon atom. The number of free backflow parameters, excluding 
cutoff lengths, is Np. The exact energy was taken from Ref. 1381 



Method Wfn. Np E (a.u.) (a.u.) % corr. en. 



HF 






-10.652399 




0.0 


Exact 






-11.055(4) 




100(2) 


VMC 


SJ 





-10.9870(3) 


0.168(1) 


83.1(9) 




BF 


214 


-11.0173(2) 


0.156(2) 


91(1) 


DMC 


SJ 





-11.0237(4) 




92(1) 




BF 


214 


-11.0348(6) 




95(1) 



TABLE VIL Slater-Jastrow and backflow results for the PP C2 molecule. The number of free backflow parameters, excluding 
cutoff lengths, is Np. The exact energy was taken from Ref. 1381 



Method Wfn. Et (a.u.) 



HF 




0.02896 


Exact 




0.233(5) 


VMC 


SJ 


0.1856(3) 




BF 


0.2051(2) 


DMC 


SJ 


0.2060(4) 




BF 


0.2093(6) 



TABLE VIII: Slater-Jastrow and backflow binding energies for the PP C2 molecule, computed using the results from Tables [Vll 
and |VII[ T he Hartree-Fock (HF) and exact energies were taken from Refs. HDll41l and 1441 and the references therein. The data 
in Ref. '38' can be used to estimate an approximate binding energy of 0.215(6) a.u.. 



4. PP diamond 

We have also studied PP carbon diamond with the experimental cubic lattice constant of 6.741 a. u.j^ representing 
the solid by a 2 x 2 x 2 supercell containing 16 atoms subject to periodic boundary conditions. Diamond is an 
insulator with a large band gap, and therefore we expect the single-determinant nodal surface to be quite accurate. 
We parametrized our backflow function using N^j = 8, = 8, and A'ee = N^.^ = 2, allowing all parameters to be spin 
and spin-pair dependent. The cutoff lengths were optimized, and they went to the maximum allowed values. The 



results in Table IX show that backflow gives a substantial reduction in the VMC energy of 0.0131(2) a.u. per atom 
[0.356(5) eV per atom], which is accompanied by a reduction in the variance by a factor of nearly two. The reduction 
in the VMC energy of diamond from using backflow is somewhat smaller than that obtained in the dimer [0.411(5) eV 
per atom], and substantially larger than that in the atom [0.147(3) eV per atom]. This may arise from the fact that 
the backflow functions in diamond arc quite long ranged and cover several atoms. Backflow reduces the DMC energy 
of diamond by 0.0035(2) a.u. [0.095(5) eV per atom] per atom, which is a little less than in the dimer [0.15(1) eV per 
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atom] and atom [0.106(3) eV per atom]. 



Method Wfn. Np E (a.u./prim. cell) (a.u.; 



DFT-PBE - 

VMC SJ 
BF 



DMC 



SJ 
BF 



-11.368208 

-11.3708(2) 
-11.3970(3) 

-11.40717(8) 
-11.4141(3) 



1.51(8) 
0.897(8) 



TABLE IX: Slater-Jastrow and backflow energies per primitive cell for PP carbon diamond using a face-centered cubic cell 
containing 16 atoms. The number of free backflow parameters, excluding cutoff lengths, is Np. 



We do not discuss the cohesive energy of the diamond crystal, as we would need to account for finite-size effects to 
be able to compare with experimental data. Within VMC, the energy gain per atom from using backflow is larger in 
the solid than in the atom, and hence the cohesive energy is substantially reduced. Within DMC, both the solid and 
the atom present a similar energy gain per atom, and the cohesive energy is not changed significantly. 



V. DISCUSSION 



A. Electron-by-electron and configuration-by-configuration algorithms 

The additional complexity of BF wave functions compared with SJ ones leads to greater computational expense 
in QMC calculations. One of the most costly operations in QMC calculations is the evaluation of the orbitals and 
their first two derivatives at points in the configuration space. The evaluation of the collective coordinates involves 
some extra cost. Furthermore, while QMC calculations with SJ wave functions require only the value, gradient, and 
Laplacian of each orbital 0, calculations with BF wave functions also require cross derivatives such as d^(j)/dxdy, 
as explained in Ref. 1181 However, the most important complicating factor arising from backfiow transformations is 
that they make each orbital in the Slater determinants depend on the coordinates of every particle. In standard 
QMC algorithms with SJ wave functions one moves the electrons sequentially in what we call the electron-by-electron 
algorithm (EBEA). Fast update algorithms are used in the EBEA to replace altered rows in the Slater determinants 
efficiently and the accept /reject step is performed on each particle separately. However, in BF calculations each 
collective coordinate depends on every electron position and therefore the fast update algorithms used in the EBEA 
are no longer appropriate, so that one must recalculate the determinants at each step using LU decomposition. 
Nevertheless, the implementation of the EBEA for backfiow wave functions can take advantage of other optimizations 
to make the algorithm more efficient, such as buffering the separate contributions to the collective coordinates, which 
we have exploited as far as possible. 

In previous fermion backffow calculation^i^IllIIIlthe electrons have all been moved together and a single accept/reject 
step has been performed, in what we call the configuration-by-configuration algorithm (CBCA). We have compared 
the efficiency of the EBEA and CBCA. The relative efficiency of the EBEA and CBCA depends on the computational 
costs of moving the electrons and the correlation time of the local energies, which is proportional to the number of 
moves of all the electrons required before the local energies are uncorrelated. Let A and B be two calculations for the 
same system, identical except for the use of different sampling algorithms. We define the relative efficiency 7 of A 
and B as 

7(A,B) = ^, (14) 

where t is the CPU time and cr is the standard error in the mean energy.!^ 7 represents the ratio of the time required 
to achieve a fixed error in the mean energy in calculation A to that required in calculation B, and is hence appropriate 
for comparing the efficiency of the two algorithms. 

In Table |X] we report results for the systems studied in this paper. For each system, the EBEA and CBCA time 
steps were chosen so that the same proportion of proposed moves were accepted: in VMC the target acceptance ratio 
was 50%, which corresponds to fairly efficient sampling, and in DMC it was around 99.5%. The correlation times 
for the CBCA are considerably longer than for the EBEA. The ratio of the correlation time of the CBCA to that 
of the EBEA (the "correlation time ratio" or CTR) appears to increase roughly linearly with the number of atoms 
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(for example, compare the AE Li atom and Li2 molecule, and the PP C atom, C2 molecule, and diamond), or with 
the number of electrons. 7(CBCA, EBEA) is larger than unity in all cases except the BF-VMC calculation of lithium 
atom, so the EBEA is generally found to be more efficient than the CBCA. 7(CBCA,EBEA) is larger for SJ wave 
functions than for BF ones because for SJ wave functions and the EBEA one uses fast update algorithms. 



System 


TVe Wfn CTR. 


7VMC 


7DMC 


HEG (r, = 1.0) 


54 


SJ 


70 


44.0 


7.5 






BF 




34.0 


1.2 


AE Li atom 


3 


SJ 


9 


4.6 


1.8 






BF 


I 


0.43 


1.9 


AE Li2 molecule 


6 


SJ 


15 


15.4 


5.0 






BF 




9.9 


3.8 


AE C atom 


6 


SJ 


9 


5.1 


7.9 






BF 


3 


2.4 


3.7 


PP C atom 


4 


SJ 


3 


3.3 


3.7 






BF 


2 


2.3 


1.3 


PP C2 molecule 


8 


SJ 


10 


6.9 


3.5 






BF 




6.9 


1.5 


PP C diamond 


64 


SJ 


80 


99.0 




(2x2x2) 




BF 




39.0 





TABLE X: Comparison of the EBEA and CBCA for SJ and BF wave functions. Key: Ne is the number of electrons; "CTR" 



is the ratio of the correlation time in the CBCA to that in the EBEA; 7 is 7(CBCA,EBEA), as defined in Eq. (14 1. Where 
a separate "CTR" for the BF wave function has not been reported, it is because it was found to equal that of the SJ wave 
function. 



Apart from the tests reported in this section, all the VMC and DMC results reported in this paper have been 
obtained using the EBEA. 



B. Computational expense of backflow calculations 

We now investigate the relative cost of BF and SJ calculations. The additional computational expense of each step 
in a BF calculation is offset by the fact that BF wave functions are generally more accurate than SJ ones, so that the 
variance of the energy is smaller, and consequently the number of statistically independent local energies required to 
achieve a given error bar in the mean energy is also smaller. 

Let A and B be two calculations for the same system, identical except for the use of different wave functions. We 
define the time ratio as r(A,B) — tx/t-e, the squared-error ratio as e(A,B) = cr^/cg, and the relative efficiency as 
7(A, B) = r(A, B)e(A, B), where t is the CPU time and a is the standard error in the mean energy. r(A, B) represents 
the relative expense per move of calculation A with respect to calculation B or, equivalently, the relative expense of 
generating a fixed number of configurations. The latter is relevant to the wave- function optimization procedure, as 
the number of configurations used should, if anything, increase with the number of parameters in the wave function. 
7(A,B) measures the relative ability of calculation A to produce a total energy to a desired degree of certainty with 
respect to B. 

In Table |XI| we compile the BF-to-SJ ratios r, e, and 7 for each calculation. For the HEG at = 20 we report 
BF3-to-SJ3 ratios instead. The performance of backflow in the HEG is impressive: backflow not only improves the 
energies but also makes the calculations less costly! 

The lithium atom is another example of improved efficiency. According to Table |TT] the SJ-DMC and BF-DMC 
energies are equal, so it does not seem advantageous to use backflow at all in this system. However, due to the 
BF-VMC energy being so close to the BF-DMC value, the variance of the BF-DMC run is enormously lowered and 
the CPU time is reduced to 25% of the time taken by the SJ-DMC iviaW 

In all cases except PP diamond, 7 is less than 3 in VMC, and 6 in DMC. However, the crystalline PP calculations 
become significantly more expensive when backflow is used. A great part of this increase is due to the computation 
of the nonlocal PP energy, which involves several evaluations of the wave function (12, in this case) for each electron 
and each ion, every time the local energy is computed. 
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System 


A^e Method 


r 






e 




7 


HEG {rs = 1.0) 


54 


VMG 


2, 


.9 


0, 


.18 


0, 


.52 






DMG 


4, 


,9 


0, 


.15 


0, 


.75 


HEG {rs = 20.0) 


54 


VMG 


1, 


.4 


0, 


.48 


0, 


.67 






DMG 


2, 


.0 


0, 


.15 


0, 


.28 


AE Li atom 


3 


VMG 


2, 


.3 


0, 


.52 


1, 


.2 






DMG 


4, 


.3 


0, 


.06 


0, 


.25 


AE Li2 molecule 


6 


VMG 


3, 


.9 





.71 


2, 


.8 






DMG 


8, 


.3 


0, 


.71 


5, 


.9 


AE C, atom 


6 


VMG 


3 


5 


0, 


.69 


2 


.4 






DMG 


5 


9 





.41 


2 


4 


PP G atom 


4 


VMG 


3 


I 





.79 


2 


4 






DMG 


2, 


.8 


0, 


.73 


2 


.1 


PP G2 molecule 


8 


VMG 


3, 


.9 


0, 


.65 


2, 


.5 






DMG 


3, 


,3 


0, 


.47 


1, 


.6 


PP G diamond 


64 


VMG 


27, 


,0 


0, 


.31 


8, 


.3 


(2x2x2) 




DMG 















TABLE XI: Data from the timing tests performed on different systems. Key: A^e is the number of electrons; t, e, and 7 are 
r(BF,SJ), e(BF,SJ), and 7(BF,SJ), as defined in the text. 



C. Backflow and nodes 

HF nodes have been compared with either exact or very accurate nodes in a number of studies. l^aiMlMlsilSS] ^ j^^s 
been found that the HF wave function often has too many nodal pockets for the ground states of atoms with four or 
more electrons. It is conceivable that coordinate transformations could modify the number of nodal pockets of a wave 
function. However, we believe this to be unlikely for the backflow transformation presented in this paper, because 
this would require the backflow displacement field to be discontinuous at very specific configurations, or exhibit other 
unusual features. The development of a general backflow transformation with the appropriate discontinuities to correct 
HF nodes, which we have not attempted, seems likely to be a tremendously difficult task. 

We now illustrate graphically how our backflow transformations changes nodal surfaces. Note that the figures 
described below are single projections of high-dimensional nodal surfaces, from which almost no useful conclusions 
regarding the full topology of the nodes can be extracted. Two-dimensional projections of the HF and BF nodes for 
a two-dimensional HEG are depicted in Fig. [8] at two different densities. The effect of backflow on the nodes is much 
more pronounced for the low-density HEG. For an unpolarized system, the nodal changes should be larger than those 
seen in Fig. |8]at all densities. Some regions of these plots suggest that the displacement of the nodes due to backflow 
is largest at points where the curvature of the nodal surface is large, away from electron-electron coalescences. There 
are a number of avoided crossings in Fig.[8](3 at = 0.5, and 6 at = 10) whose connectivity (in the projection) is 
modified by backflow. 

Three-dimensional projections of the HF and BF nodes of the AE carbon atom are compared in Fig. |9] The nodes 
are substantially modifled by the introduction of backflow. New nodal regions appear in this projection because the 
electron being moved "pushes" the other electrons (via the backflow transformation) through the nodal surface of the 
HF wave function. 



VI. CONCLUSIONS 



We have devised an inhomogeneous backflow transformation for systems consisting of electrons and either nuclei 
or ions represented by pseudopotentials. We have applied our backflow transformation to single-determinant Slater- 
Jastrow wave functions for the HEG and for atomic, molecular, and solid systems. In each case backflow gives a 
substantial reduction in the VMC energy, and a smaller reduction in the DMG energy. 

The homogeneous backflow transformation reduces the variance of the VMG energy of the HEG by a factor of 
about 4, which is the largest such factor we have encountered, and we believe that our backflow wave functions for the 
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(b) Nodes in a 2D HEG at = 10. 

(a) Nodes in a 2D HEG at = 0.5. 



FIG. 8: (Color online) Nodes encountered when moving one of the electrons of a two-dimensional HEG of 101 like-spin electrons 
at two different densities. The HF and BP nodes are in black and orange, respectively. The green circles indicate the positions 
of the remaining 100 electrons, at which the nodes are required to remain fixed. The backflow wave functions were obtained by 
variance minimization; the energy reductions from SJ-VMC to BF-VMC at rs = 0.5 and 10 were 0.0007(1) a.u. /electron and 
0.00005(1) a.u. /electron, respectively. 

HEG are very accurate. VMC retrieves more than 99.5% of the DMC correlation energy in the density range studied 
(rs = 0.5-20). The effects of backflow on the nodes increase with r^, even though the additional percentage of the 
correlation energy retrieved in VMC decreases with r^, implying that the energies of dilute HEGs are less sensitive 
to the nodal structure of the trial wave function than those of denser systems. 

Although backflow works very well in the HEG, as previous studies have already concluded, we find that purely 
homogeneous backflow transformations give poor results when atoms are present, as we demonstrated for the AE 
lithium dimer and the AE carbon atom. However, in these cases inhomogeneous backflow transformations can 
improve the wave functions substantially. 

For the AE lithium atom the HF nodal surface of the SJ wave function is essentially exact. Although in this 
case backflow cannot improve the DMC energy, it gives a very accurate VMC energy. This shows that backflow 
transformations can improve the wave function away from the nodes as well as improving the nodal surface itself. The 
quality of the SJ and BF wave functions for the AE lithium dimer are much lower than for the atom, and consequently 
the binding energy of the dimer is underestimated. The wave function and nodal surface of the AE lithium dimer 
can be substantially improved by using several determinants, "^"^ but it appears that only modest improvements can be 
obtained using backflow. 

Backflow reduces the VMC energy of the AE carbon atom by about 49% of the correlation energy missing at the 
SJ-VMC level, but at the DMC level the improvement is smaller; the BF-DMC energy is only 18% closer to the 
exact value than SJ-DMC. Backflow makes a more significant improvement to the DMC energy of a PP carbon atom 
than the AE carbon atom. The PP and AE carbon atoms are also cases where substantial improvements to the 
wave functions can be obtained by using several determinants. This indicates that the SJ nodal surfaces of these two 
systems need a more drastic correction than backflow transformations can provide. 

When the initial nodal surface is reasonably accurate, backflow does an excellent job in improving the VMC energy 
and correcting the remaining errors in the nodal surface, as was seen in our study of the HEG and AE lithium. 
However, when the initial nodal surface is intrinsically poor, as is the case, for example, with the HF nodal surfaces 
of the carbon atom and dimer, backflow is apparently incapable of making the gross changes to the nodal surface 
required to correct the flaws, although it normally lowers the VMC and DMC energies somewhat. We do not believe 
that our backflow transformation is capable of changing the number of nodal pockets of the starting wave function. 

The cost of using BF wave functions can be substantial, but we have given evidence that the expense relative to that 
of using SJ wave functions increases smoothly with the number of atoms in the system. Backflow transformations, like 
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(a) HF nodes of a carbon atom. (b) BF nodes of a carbon atom. 



FIG. 9: (Color online) HF and BF nodes encountered when moving one of the (majority spin) up-spin electrons of an AE carbon 
atom. The blue circle corresponds to the position of the nucleus, the red upward-pointing triangles indicate the positions of 
the remaining up-spin electrons and the green downward-pointing triangles indicate the positions of the down-spin ones. The 
HF node consists of a (seemingly) infinite sheet with a bubble attached to it, which contains the nucleus. Backflow slightly 
modifies this node, and adds three large lobes (detached from one another; all intersect the HF node) and a small bubble next 
to a down-spin electron. 

Jastrow factors and unlike niultideterminant expansions, are compact parametrizations, meaning that the number of 
parameters required to retrieve a given fraction of the correlation energy increases only slowly with system size. This 
can be seen by comparing the number of backflow parameters that we have used and the energies we have obtained for 
PP carbon atom, dimer, and diamond. We have found that it is much more efficient to move electrons one at a time 
(the EBEA) than to move all the electrons at once (the CBCA) , as has been done in previous backflow calculations. 
The reason for this is that the correlation time of the energy is considerably shorter with the EBEA. It is important 
to use the EBEA for large systems, as the CBCA-to-EBEA ratio of correlation times seems to increase linearly with 
the number of electrons. 

BF-VMC energies are normally significantly lower than SJ-VMC ones, and therefore BF-VMC might be a useful 
alternative to a (normally more expensive) SJ-DMC calculation. The use of more accurate trial wave functions 
improves the statistical efficiency of VMC and DMC calculations. The variance of the local energies encountered 
in a DMC calculation is approximately proportional to the error in the VMC energy, and when backflow leads to a 
significant reduction in the VMC energy it also improves the statistical efficiency of DMC calculations, even when 
backflow improves the DMC energy only slightly. The improved trial wave functions could also be useful in DMC 
calculations of quantities other than the energy, which are normally more difficult to obtain accurately than the 
energy. 

Backflow would appear to give significant improvements in trial wave functions for a wide variety of systems, in- 
cluding various different atoms, and small and large systems. In the present work, we have applied the inhomogeneous 
backflow transformation to single-determinant Slater-Jastrow wave functions only, but it can be combined with multi- 
determinant wave functions, and we will report on such calculations elsewhere It can also be combined with pairing 
wave functions.'^ We believe that inhomogeneous backflow transformations will play an important role in improving 
trial wave functions for use in VMC and DMC calculations. 
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APPENDIX A: CONSTRAINTS ON THE BACKFLOW PARAMETERS 



1. Cusp conditions 

The Kato cusp conditionpEs! (KCC) are enforced so that the local energy is finite when two electrons or an 
electron and a nucleus are coincident. For SJ wave functions it is common practice to impose the electron-electron 
KCC (EKCC) by constraining the parameters in the Jastrow function, and the electron-nucleus KCC (NKCC) by 
constraining the orbitals in the Slater determinant. The backflow transformation can alter the nature of the cusps, 
but we have chosen to constrain the backflow parameters so that they do not modify the KCC as applied to the 
Slater-Jastrow wave function.-^ 

Let i and j be two different electrons in the system. To satisfy the EKCC, we require that the total backflow 
displacement ^j, has a well-defined gradient (i.e. it should be cuspless) when Vij ^ if i and j are distinguishable 
particles, and have zero gradient when if i and j are indistinguishable. Thus, the e-e term is affected by these 

constraints only if i and j are like-spin electrons, in which case the EKCC are satisfied if LfjCi = Ccq. 

Let / be a nucleus in the system. To satisfy the NKCC, we require that the total backfiow displacement, has a 
well-defined gradient when ru 0, and that it is zero when rn — > if / is an AE atom. The NKCC are satisfied if 
L^jdi j — Cdoj for all /, and in addition, cIqj = 0, if / is an AE atom. 

The constraints on the e-e-n functions, some of which only apply to those functions centered on AE atoms, are as 
follows. (We omit the / index in the parameters for clarity.) 

• There are 3 (A'co + N^n + 1) constraints from the NKCC, 

(l-\-m—a) (/c+m— a) (k+m—a) 

{CcpQim — L$(/3i;„i) = ^ {CcpkOm " Lq,ipklm) — ^ (C'^feOm " Lq,9klm) =0 V a . (Al) 
Z,m fe,m k,m 

• There are 27Von + 1 constraints from the EKCC, 

(k+l=a) 

dkii^O Va, (A2) 



k,l 

and 2Nen + 1 extra constraints for like-spin electron pairs. 



{k+l=a) 

Vkii = Va. (A3) 



k,l 



[AE only] There are 4 {Ncc + Ncn) + 2 constraints on (pkim, 

(l+rn—a) (l+m—a) (fc+m— a) (fc+m— a) 

"^Oim = Y '^'POlni = Y '^'=0™ = X! "^VfeOm = V Of . (A4) 
/,m l,m k.m k.m 



[AE only] There are 3 (A^cc + -^cn) + 2 constraints on 9^ 



klr. 



a) (/+7n— a) (fc+m— a) 

Y ^0'™ ^ X! ™^Oim = Y "^^feom = V a . (A5) 
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These constraints form an indeterminate system of homogeneous algebraic linear equations for the e-e-n parameters. 
Hence, a subset of the parameters can be put in terms of the rest. This subset can be determined from the "free" 
parameters by putting the constraints in matrix form and using Gaussian elimination. This procedure is the one 
described in Ref. where it is applied to the parameters in the e-e-n term of the Jastrow factor. 



2. Constraints for irrotational backflow 



In the derivation of homogeneous backflow in Ref. [2T]it was suggested that the backflow displacement should satisfy 
= V iY , where Y = ^(R) is an object called the backflow potential. This equation is already satisfied by both 

the e-e and e-n terms, by definition, and it can be imposed on the e-e-n functions by using an appropriate set of 

constraints. From x = 0, it follows that 



d 
' dru 



d 



for all i, j, and /, and all r, 

(C + k) ipk,i,, 
while for C = 0, 



ij, rn, and rjj. For C > 0, this results in the equation 

-1 - i$ (fc + 1) ¥?fc+l,/,„i_l - (to + 1) 9k-2,l,m+l + (to + 1) 9k-\,Lrn+l = 
(fc -I- 1) (Pk+l.l,,n-l - (to -I- 1) 6k-l,Lm+l = . 



(A6) 



(A7) 



(A8) 



In both cases, < fc < iVcn + 2, < Z < A'en, and < m < iVee -I- 1, and parameters with indices out of the allowed 
range are to be taken as equal to zero. The I index has been omitted for clarity. 

The application of these constraints results in a reduction in the number of free parameters by more than one half, 
as one would expect, because an equivalent backflow displacement would be obtained by parameterizing the scalar 
field Y and computing its gradient, whereas we use two scalar fields in the full e-e-n term. 



APPENDIX B: ZEROING THE BACKFLOW DISPLACEMENT AT AE ATOMS 



When AE atoms are present, the NKCC cannot be fulfilled unless the backfiow displacement at the nuclear position 
is zero. This can be obtained by applying smooth cutoffs around such atoms. In this scheme, an artificial multiplicative 
cutoff function girn) is applied to all contributions to the backflow displacement of particle i that do not depend on 
the distance rn to the AE atom /. This includes the homogeneous backflow displacement and the inhomogeneous 
contributions centered on each atom J ^ I. 

The g{rii) function must go to zero at rn — > and become unity when ru is equal to or greater than a threshold 
Lgj. For the local energy to be well-defined, we require that g{rii) and its first two derivatives be continuous at 
ru — Lgj, and to fulfill the NKCC correctly, g{rii) and its first derivative must go to zero at ru = 0. The simplest 
g{rii) obeying these conditions is the fourth-order polynomial. 



Lg.I 



6- 



L 



(Bl) 



which we have used in our calculations. Although it is perfectly possible to optimize the Lgj, we have used fixed 
values for simplicity: 1 a.u. in the AE atoms and about half the interatomic distance in the AE Li2 molecule. 
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The C-th derivative of the wave function will be discontinuous at r = L. In particular, its Laplacian, used in the computation 
of the kinetic energy, is discontinuous at r = L if C < 3. In this work we have only considered C — 2 and C = 3. 
H is the Hamiltonian operator of the system, and the total local energy is defined as El = 'if^^H'^/r- 

The VMC variances for the HEG at = 1 reported in Table II of Kwon et a?!'-® have been confirmed by the authors to be in 
error; the true values are a factor of 10 smaller. These data were later used in Table II of Holzmann et aZ.,'^ who corrected 
the mistakes, except for the variance of the SJ calculation. We have compared our data with the corrected values. 
Two configurations are in the same nodal pocket if there exists a continuous path between the two along which the wave 
function does not change sign and is not equal to zero. Nodal pockets are bounded by nodal surfaces, which determine the 
shape and number of the former. „„„ 
The variance of the local energies encountered during a DMC calculation is approximately proportional to -Evmc — ^■dmc-'^^'^ 
In principle, it would be possible to apply the KCC to the Jastrow and backflow parameters together. However, the resulting 
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constraints are configuration-dependent and involve orbital derivatives, making this approach difficult. 



